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Q ' Abstract 

We study power corrections to the differential thrust, heavy mass and related event shape 
distributions in e^e~ annihilation, whose values, e, are proportional to jet masses in the two- 
^ ' jet limit, e — > 0. The factorization properties of these differential distributions imply that 

j^ ' they may be written as convolutions of nonperturbative "shape" functions, describing the 

emission of soft quanta by the jets, and resummed perturbative cross sections. The infrared 
shape functions are different for different event shapes, and depend on a factorization scale, 
but are independent of the center-of-mass energy Q. They organize all power corrections 
of the form l/(eQ)", for arbitrary n, and carry information on a class of universal matrix 
elements of the energy- momentum tensor in QCD, directly related to the energy-energy 
correlations. 
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1 Introduction 

Infrared safe event shapes in e"'"e~ annihilation are finite order-by-order in perturbation theory 
and can be calculated at large center-of-mass energy Q, in powers of as{Q). Nonperturbative 
corrections decrease as powers of Q relative to the perturbative series. In most cases such 
corrections are readily observed in the data up to the highest available energies. They can be 
parameterized as Xp/Q^, with scales Xp and exponents p depending on the shape in question, and 
they must be taken into account in precise measurements of as{Q) based on event shapes |T|, 0. 

Some time ago, it was shown that perturbation theory itself suggests the exponents p (but 
not the scales Xp) for the leading power corrections to the mean values of various event shapes 
[0, ^, |, |[. The simplicity of these results, and the successes of fits to the data based upon them, 
were somewhat surprising, because previously nonperturbative effects were estimated primarily 
by comparison to event generators. 

Although perturbation theory cannot predict the magnitudes of the scales Xp, it is natural to 
entertain the possibility that the leading power corrections to event shapes might be controlled by 
a set of nonperturbative but universal parameters 0, |^, ^ . The latter have been parameterized 
as moments of a universal effective coupling over low momentum scales [^ , and further analyses 
incorporate effects of the effective coupling beyond lowest order 0. It is not clear, however, 
that long-distance effects can be expressed fully in terms of low powers of the running coupling, 
extrapolated to scales of order Aqcd- We therefore believe that a more general viewpoint will be 
desirable, to further reveal the underlying QCD dynamics. 

In this paper, we suggest a formalism for power corrections to IR safe event shapes, based on 
the factorization of short-distance perturbative and long-distance nonperturbative effects. The 
factorization formalism can be thought of as a generalization of the operator product expansion 
(OPE), applicable to weighted cross sections. For event shapes, we encounter a new set of 
nonperturbative distributions of soft radiation, which we shall term "infrared shape functions" 
[0. As we observe below, our infrared shape functions generalize the inclusive (light-cone) 
distributions that have been introduced to describe semileptonic b quark decays near the edge of 
phase space [0. They admit operator definitions, and are related to correlation functions of the 



energy-momentum tensor, which have been discussed in a number of related contexts ||12| , |T^, [14 



The moments of the infrared shape functions define the nonperturbative scales Xp that determine 
the mean values of various event shapes. Previous descriptions of power corrections reemerge as 
particular choices for these functions. 

Rather than analyzing only mean values, we will treat differential event shape distributions. 
The structure of power corrections in this case is more complex, and more interesting. For 
example, it has been shown for the thrust that when t = 1 — T is small, but substantially 
larger than Aqqd/Q, the 1/Q power corrections produce a shift in the resummed perturbative 
distribution 0,0,0, 

^ + 0(l/(tQ)2), (1) 



dt dt 

with Ai the same scale that parameterizes the 1/Q correction to the mean thrust, (t) ~ Xi/Q. 
This approach improves the fit of the resummed perturbative formula for low values of t and 
related event shapes. At the same time, Eq. (|I|) must break down for very small values of 
t ~ Aqcd/<5, corresponding to small invariant mass jets in the final state, because in this 
region all power corrections of the form l/(t(5)", with a = 1,2, .. . become equally important. 
We shall observe that the factorization properties of event shape distributions enable us to 



organize all powers of l/{tQ), by expressing the distributions as convolutions of nonperturbative 
infrared shape functions with resummed perturbative spectra. The scale for this factorization, 
/i ^ AqcD) will act as an infrared (IR) cutoff in the perturbative expressions. Although the 
physical distributions are independent of /i, the shape function and the resummed spectrum for 
each event shape depend, as usual, on the factorization scale. The resulting expressions for the 
differential distributions are applicable over a wide interval of the event shapes, and for moderate 
t they take the simple form of Eq. (0). 

We shall discuss two representative event shapes in some detail, the thrust, t, and the heavy 
jet mass distribution, p. These two shapes are ideal testing-grounds for our proposals on the 
interplay of perturbative and nonperturbative effects in QCD, since we have fairly good control 
over their perturbative expansions, while at the same time the data indicate the presence of 
significant power corrections. 

We begin, in the following section, by showing that the kinematics of e"'"e~ annihilation near 
the two-jet limit suggests that the heavy jet mass and thrust distributions are convolutions of 
nonperturbative and perturbative components. In Sec. 3 we show that the same convolution form 
emerges from known factorization properties of the resummed cross sections in the two-jet limit. 
We go on in Sec. 4 to show how the infrared shape functions are related to the expectation values 
of products of Wilson lines, independent of the center-of-mass energy. We then show how the 
truely universal nonperturbative information in event shapes is encoded in a set of correlation 
functions of an energy flow operator, related to energy-energy correlations. 

2 The Radiation Function and Factorization 

Consider an event shape e (= 1— T , p , . . .), which takes on a value e(A^) for each e"'"e~ annihilation 
final state |A^). We represent the differential distribution da/de schematically by 

1 da{Q) ^^|^^|j(o)|0)r5(e-e(Ar))^(5(e-e(Ar))), (2) 



(Ttot de o"tot jv 

where J{x) represents the appropriately normalized electroweak current, and the sum over states 
N is taken at fixed momentum g, q^ = Q^. The second equality in Eq. (§) defines a convenient 
notation for the weighted distribution. 

The weights e{N) for the thrust and jet mass are given in the end-point region, e ^ 0, in 
terms of the total invariant masses M|. ^ flowing into right and left hemispheres defined by the 
plane orthogonal to the thrust axis, as 



t{N) 



Ml + Ml 



Q' 



/Vf2 11^2 

p{N) = -^e{Ml-Ml) + -^e{Ml-Ml). (3) 

Each soft particle with momentum £fc = {Ek,ik) in the final state contributes to these weights 
additively. Its contribution depends on whether it goes into the right {cos 9k > 0) or left (cos6'fc < 
0) hemisphere, with 6k the angle between the 3-momentum of the particle, ik, and the thrust 
axis. Neglecting terms quadratic in the energy of the soft particles, we find 

MIl = QI: i5;,(l-|cos^,|), (4) 

fc=l 



where Nh^l denotes the number of soft particles in the corresponding hemisphere. 

We now introduce a "radiation function", -R(e), in the notation of Eq. (|^), through [|T6|1 

R{e) = {6 {e - e{N))) ^ de' — ^ = 1 - d^' — fj ■ (5) 

Jo (Jtot de' Je (Ttot de' 

Here, Cmax is the maximum value of the event shape, e, and atot is the integral of the differential 
distribution over its full range, cxtot = Jo^'"'"' deda{e)/de. The function R{e) is unity at e = Cmax 
and vanishes as e — ;> 0. We expect perturbation theory to be accurate for R{e) at "large" 
e ~ Cmax; ^ud to fail when e = 0{A/Q). In the notation of Eq. @, the radiation function for 
the heavy jet mass is 



RH{p) = {9ip-piN))) = (9 p--^\9 p--^\), (6) 




with Mr and M^ given in Eq. (^). Similarly, the radiation function for the thrust is 

Rrit) = {6 it- m)) = (e[t- ^^^^) ) . (7) 

The functions Rh and Rt receive both perturbative and nonperturbative contributions. 

For such event shapes, e = corresponds to a final state with two infinitely narrow jets, and 
the resummed perturbative differential cross section vanishes in this limit, while its expansion in 



as{Q) is singular order by order |16|. For small e, the final state consists of two narrow quark 
jets of invariant mass ~ Qe^^"^ and a "cloud" of soft particles of the total energy ~ Qe. We 
expect nonperturbative corrections to occur with inverse powers of both scales. However, as e 
decreases, the energy of the soft radiation reaches a nonperturbative scale before the jet invariant 
masses, because Qe^^^ ^ Qe. Therefore, for e^^^ ^ ^qcd/Q ~ e, we may restrict ourselves to 
power corrections in the smaller scale Qe. We shall see below that it is possible to organize 
systematically power corrections of the form l/(Qe)", with a = 1,2.... With this in mind, we 
consider the limit e — *^ with Qe = fixed, and neglect corrections suppressed by additional 
inverse powers of Q. 

In this limit, the jet invariant masses M'^^ are sensitive only to the light-cone components of 
soft particle momenta in the final state, £^ = Ek{l ± cos9k)/V^, defined relative to the thrust 
axis. Assuming that nonperturbative effects are small, we may separate the contributions of soft 
particle emission to Ml and Mr as 

\rsoft 

MI^MI + srQ, MI^MI + slQ, £r,l= ^ Efe(l-lcos^fcl), (8) 

fc=i 

where from now on, M'j^^^^ refer to the perturbative contributions to the hemisphere masses 
squared, and Sr^lQ to the soft-parton, nonperturbative contributions. Because, as we now argue, 
soft particle emission may be factorized from the jets, El^r are sums over the light-cone momentum 
components of Nf^^^ soft particles, emitted into the respective hemispheres. 

The basis of factorization is that the emission of soft gluons occurs over time scales that 
are different from those involved in the evolution of narrow jets. These processes are quantum- 
mechanically incoherent, and the soft gluon distribution emitted by a pair of narrow jets at wide 



angles depends only on the direction and total color charge of the jets. (See Sec. 4 below.) 
Then, to all orders in perturbation theory, inclusive cross sections for two narrow jets in e"'"e~ 
annihilation can be written as products of separate functions for the jets and for the soft radiation, 
convoluted in the light-cone components of the soft radiation. Corrections are suppressed by 
powers of M^^/Q^. Details of the necessary reasoning are reviewed in Ref. |T^. Here, we apply 



this factorization to the limit e — ;> 0, Qe fixed, identified above. In this limit, the jets are narrow, 
but their invariant masses are still large enough for perturbation theory to be valid, while the 
soft-gluon function becomes nonperturbative. 

For the heavy jet mass, we use Eqs. (|^) and (|^) to identify the appropriate factorized expres- 
sion as 



Rh{p) = j denj deL f{€R,eL) (o ip - 



Ml ^.\J^_Ml.-^^^ , (9) 





where the subscript PT denotes the average with respect to the perturbative spectrum in Eq. (H), 
in a manner we will specify below. The factor f{eR,eL) is the nonperturbative infrared shape 
function referred to above. It represents the probability density for the total soft gluon light 
cone momentum components, Er and El, in each hemisphere. The factorization of soft dynamics 
implies that /{er, El) does not depend on the hard scale Q, up to corrections suppressed by l/Q, 
although it does depend on the cut-off /x that sets the maximal energy of particles described by 
the shape function. The /i— dependence of the shape function is compensated in (^) by that of 
the perturbative contribution. (For now, we suppress p as an argument for simplicity.) 
For the thrust, we have, with the same shape function f{ER,EL), 



Rrit) = l^^dERJ^^dEL fiSR.EL) 1 9 It 



^^-^)> . ao, 



The shape functions in Eqs. (|^) and (p!0| ) are identical precisely because of the factorization of soft 
gluon emission from the jets. Note that f{ER,Ei,) is a symmetric function of the soft radiation 
variables Er^l separately. Because Rt depends only on the sum of soft momenta, we may simplify 
it by introducing the function 

fj,[e) = dER dEL /{er, El) 5{e -Er- el) , (11) 

Jo Jo 



in terms of which 



Thus, the leading power corrections to the thrust give rise to an e/Q— shift of the perturbative 
radiation function, just as in Eq. (|l|), but now averaged with the shape function fri^)- In 
contrast, the jet mass function, Eq. (|[), does not have this property, unless 

fiSL, Er) = f\ER) f{EL) + Af{EL, Er) , (13) 

with A/ small, corresponding to independent evolution for the two hemispheres. This is not 
obvious even perturbatively, however, and Af{EL,ER) is nonzero, for example, due to the decay 
of an off-shell gluon into a pair of gluons, each moving into different hemispheres Jl 



Even though Eq. (|l^) is not a general property, the perturbative expectation values do factor- 
ize into separate functions for the jets to resummed next-to-leading logarithmic approximation 

(NIL) |T6l, 
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(14) 



We will introduce the infrared cutoff into these perturbative functions. Following Eq. (jT^), we 
for the heavy jet mass distribution, 



rpQ rpQ ( ^R \ ( ^L ^ 

Rh{p) = ] dsRJ deLf{eR,eL,^)RY \p- —,^,\ Ry \p- —,^i 



(15) 



where Ry^{p,fi) is the cutoff-dependent perturbative radiation function for a single jet with 
invariant mass pQ"^, 



RTip.i^) 




d{p) 



1-1 dp' -^ 

ma.x{p,fi/Q) (Jtot "P 



(16) 



In this expression, da^^ /dp' is the inclusive perturbative cross section for the creation of a single 
jet of invariant mass pQ^, and we have exhibited the p dependence for clarity. Compared with Eq. 
(^, an infrared cutoff is imposed on p, by demanding that pQ^ > pQ, which corresponds to the 
mass of a jet that includes radiation with energy of order p at large angles to the jet axis. In this 
way, we introduce p as the factorization scale that separates perturbative and nonperturbative 
dynamics. The perturbative radiation function, -Rj^, is "frozen" for < p < p/Q. This is 
equivalent to cutting off the perturbative spectrum at the value p = p/Q. All dynamics below 
p = p/Q is to be incorporated into the soft shape function. Other choices for the transition 
between perturbative and nonperturbative regions are possible, but this is relatively simple to 
implement numerically |l^. Note that -R"(0,p) = RY{p/Q, p) includes the contributions only 
of virtual soft gluons with energies above p. Therefore, because p is nonzero, infrared divergences 
do not appear and Ry^{0,p) does not vanish, although it is small. 
For the thrust, the corresponding expression for -Rr(^) is 



ftQ 

Rrit) = / de Me) R 
Jo 




(17) 



where Rj^(t',p) is the perturbative radiation function corresponding to independent evolution 
of two jets, the sum of whose invariant masses is bounded by t'Q^: 
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Here, the IR cut-off is implemented on the invariant mass of each jet and on their sum. In the 
remainder of this section, for simplicity we shall again suppress the argument p in the cut off 
perturbative and infrared shape functions, although they are always present. 

The thrust and heavy jet mass distributions themselves are found by taking derivatives of their 
corresponding radiation functions, Eqs. (|l^) and (|l^), respectively. Consider first the thrust, for 
which 



1 da dRT 



dt 



o"tot dt 



ftQ 

QfT{tQ)RV{0)+ / deMe) 
Jo 



1 daPT^ it 



crUt 



dt 



(19) 



where 

crfot dt dt 

Here, the first term corresponds to freezing perturbative real soft gluon radiation, and replacing 
it with the nonperturbative shape function fxitQ) (multiphed by the cut-off Sudakov form factor 
i?j?'(0)). The second term is the perturbative spectrum smeared with the shape function. The 
resulting form, Eq. (|19D, has a simple physical meaning. The leading l/(tQ)— nonperturbative 
effects act independently of perturbative branching, and the physical distribution is obtained by 
convoluting the perturbative spectrum da^'^ /dt with the nonperturbative infrared shape function 
fri^)- Although da^'^ /dt and /t(^) depend separately on the factorization scale /i, the physical 
cross-section da/dt is /i— independent. As emphasized above, the nonperturbative shape function 
/r(e) depends on /i but not on the hard scale Q. It takes its maximum value for e ~ /i and 
rapidly vanishes for larger e. In contrast, the perturbative spectrum starts at t ~ fi/Q and 
extends to tmax- For very small < t < fi/Q the shape of the thrust distribution is governed 
entirely by fritQ)- 

For the heavy mass distribution we differentiate Eq. (^) and use the symmetry of the shape 
function, /(£L,£i?) = f{£R,£L) to obtain 

1 da [pQ 1 day (p - f^ 

— -r = QfHipQ,pQ)RV{0) + / defH{e,pQ) ——^ ^ 

CJtot dp Jo crfot dp 

where we have introduced the function J'h, 




Me, pQ)=2 de' f{e, e')Ry{p--] . (22) 



Although the shape functions /t and fn in Eqs. (|T9D and ( pTj) are related to the same distribution 
f{£L,£ii), in general they are different. In contrast with fx, the shape function fn depends both 
on the center-of-mass energy and the shape variable. For pQ > p we may expand the shape 



function (^) in powers of l/{pQ) as 

1 dlnRYip) r,,,,,^ . , ,^( 1 



fH{e.pQ)=2RY{p) 



Jo pQ dmp Jo \{pQ) 



2 



(23) 



Comparing Eqs. (|T^) and (pT]), we see that the leading power corrections to the thrust and the 
heavy mass distribution cannot be described by a simple shift, associated, for instance, with a 
universal effective coupling constant. It is, however, relatively straightforward to use Eqs. (|T^) 
and (^) phenomenologically [^ . 



3 Soft Radiation in Perturbation Theory 

Let us now show how the principal features of nonperturbative shape functions emerge in re- 
summed perturbation theory. By deducing the Q-dependence of nonperturbative contributions 
from ambiguities in the perturbative expansion 0, ^, |, ^ |^, ^, we will show that leading non- 
perturbative contributions have the properties of the infrared shape functions identified in Sec. 
2. Specifically, they enter the cross section in convolution with perturbative cross sections, and 
they are independent of Q. To be specific, we will concentrate on the thrust distribution da/dt. 



We start by recalling that the cross section da/dt computed in the end-point region t ~ by 
resumming an infinite number of soft gluon emissions, exponentiates under a Laplace transform, 

rfte-'^*^ = e-^(^'«). (24) 



(Ttot "'0 dt 



This property has been derived in a variety of ways |T6|, |T9|, ^ , and is related, in particular, to 
the independent evolution of the opposite-moving jets. The exponent S^u, Q) in Eq. (p^ is of 
the general form 



S{y,Q) = f- (l-e— 
JO a ^ 



aQ2 77,2 



= SpT{iy,Q,ij) + Snp{u,Q,jj), (25) 

where in the absence of an IR cut-off the integration over soft gluon momenta extends to ar- 
bitrarily small values of k±. As indicated in the second equality, we will separate perturbative 
and nonperturbative parts of S by introducing a factorization scale fi. Heuristically, the contri- 
butions of gluons with k± > fi are to be absorbed into Sp^, and those with k± < fi into a new, 
nonperturbative exponent S'np- 

Consider first the perturbative exponent Sp^{v,Q, fi). Its properties depend on the value of 
the Laplace transform parameter z/, which in turn is conjugate to the thrust variable, 1/z/ ~ t. 
For 1/z/ > fi/Q one can replace (1 — e~'^") — > 9{a — l/z/e~'^^) in (E3) and derive an expression 



for Sp^ that resums all leading and next-to-leading logarithms of z/ |16|. 

For a > 1/z/, the coupling in the term B{as{aQ^)) is always at a perturbative scale, and S'np 
comes entirely from the F term in Eq. (p5|). In general, S'np depends on two scales, u/Q and 
jy/Q^. Formally expanding in powers of z//Q, we find 
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where we have neglected all powers of z//Q^, which are suppressed by powers of 1/Q relative to 
those we have kept. In the final line of Eq. (p6|), we absorb the ambiguous integrals over soft 
gluon momenta into the definitions of nonperturbative but Q-independent parameters A„(/i). 
Since the variable v is conjugate to t, the final line of ( pUf ) organizes all power corrections in 
l/{tQ). This infinite set of nonperturbative parameters A„(/i), defines a single, Q-independent 
function through 

/■oo 

/ de exp(-z/£:/(5)/T(£:;Ai) =exp(-SNp(Q/z^,/i)) , (27) 

Jo 

where, as the notation suggests, we may identify this function with the infrared shape function 

for the thrust cross section. 

Comparing the expansions of the two sides of this relation in powers of v /Q we derive the 

sum rules 

I de frie; fi) = 1 , I de e frie; fi) = \i , I de e^ frie; fi) = \l - \2 , ■ ■ ■ (28) 



Eqs. (^) and (|2q) give the Laplace transform of the thrust distribution ( plf ) as a product of 
perturbative and nonperturbative functions |3 Replacing the nonperturbative exponent in Eq. 
( p5D by its expression ( pT]) and performing the inverse Laplace transformation, we obtain the 
thrust distribution da/dt in the form of a convolution of the perturbative Sudakov spectrum 
da^^/dt and a Q-independent nonperturbative shape function fT{e,fi), as in Eq. (p!9|). 

We emphasize that this formalism is designed to treat t > /i/Q, or equivalently u < Q/^. 
Beyond this region, the soft function acquires Q-dependence through terms that are neglected 
here, including those associated with B in Eq. (^). Similarly, for u > Q/fi, it is necessary to 
go beyond the NLL approximation in S'p^. We shall defer the treatment of these extentions to 
future work. 

Summarizing our analysis of soft gluon resummation, we conclude that, in accordance with 
our expectations, the factorized expression for the thrust distribution, Eq. (|1^), resums large 
perturbative logarithms and organizes all power corrections of the form [l/tQ)^. The latter are 
described by a single infrared shape function /r(£) that takes into account the contribution to the 
thrust of soft gluons with energy smaller than /i. It is straightforward to repeat similar analysis 
for the heavy mass distribution to arrive at Eq. (^Tj). 

4 Correlation Functions 

We now turn to the field-theoretic content of the shape functions. Because the IR dynamics 
of the soft gluon radiation described by /{erjEl) and fri^) factorizes from the quark jets, we 
can apply the eikonal approximation, and replace the jets by eikonal lines in the computation 
of the soft emission. Let p+ and p_ be the light cone directions defined by the quark jets. The 
corresponding eikonal, or "Wilson", lines, which serve as sources of soft radiation, may be written 
as $±(0) = P exp [ig J^ ds Pj.Afj_{p± ■ s)]. In this approximation, the shape distributions become 



1 daiRJe) 
(Jtot de 



Y: \{N\WmO)\'5{e-e{N)), (29) 
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where W{0) = $+(0) ($_(0)) is the total eikonal phase of the two quark jets 0, and where the 
sum is over final states |A^). 

As it stands, Eq. (|29| ) requires renormalization, which we can carry out with scale fi serving as 
a UV cut-off on soft gluon momenta contributing to S-^p. Recalling that /i appears as an IR cutoff 
in the perturbative expansion of the infrared factor in Eq. (^^, we may identify daiB^{e)/dt in 
( p9D directly with the infrared shape function for observable e. For example, the shape function 
f{eL,£R) corresponds to the particular choice of e as the pair of light cone variables {eL,£R) 
defined in Eq. (||). 

To express the shape function as a vacuum expectation value, we introduce an operator S{n) 



that measures the density of energy flow in the direction of unit vector n [^. The operator £{n) 
is defined by its action on the final states of A^ soft gluons with the energy Ei and momentum 
Eifii 

N 

£{n)\N) = J2Sicos 9 - cos 9,) 5{ip - ipi) Ei\N) , (30) 



^Although this result was obtained without use of the OPE the same property holds for quantities which admit 
the OPE (for example, DIS structure functions for x ^ 1). 



with 9 and cj) the spherical angles defining the position of unit vector n with respect to the 
thrust axis. It admits a representation in terms of the energy-momentum tensor at infinity 
0, |13|, |I^. Using the definition (pO]), we find the following operator expression for the shape 
function f{eL,en}, 

f{eL,en) = {0\W\0)6 (^Bl - J dnwL{n) ^(n)) 6 (^Br ~ J dnwR{n) ^(n)) H^(0)|0) , (31) 

and, as a result of Eq. ([TT|) , for fri^), 

Mb) = {0\W\0)6 (b - f dnwrin) S{n)] W{0)\0) . (32) 

In these expressions, the integration is over unit 3-vector dn = d^n6{l — \n\) = d{cosd)dip, and 
the weight functions, wt = wr + wl, are given by0 

'W^i?,L(^) = (1 — I cos^l) 6(±cos^) , wrin) = 1 — \cos9\ . (33) 

The factorization scale fi has the meaning of the normalization point of operators entering (|32|) . 
We notice that the weight functions suppress the emissions of gluons coUinear to quark jets, 
w{n) ^ as cos 6 -^ ±1, and that therefore the shape functions do not suffer from coUinear 
singularities. 

The function fri^) cannot be considered as a universal distribution for the final states of e"''e~- 
annihilation, even in the two-jet limit, since it reflects the choice of the weight wt- However, 
taking its integer moments with respect to b, 

{b^) = dBB^ friB, /i) = / n ^^i wririi) g{ni . . . n^v) , (34) 

we encounter the Green functions Q given by multiple correlators of the energy flow operators 
£(n), 

g{ni ...nN;f^) = {0\W\0)S{ni) . . . ^(n^)iy(0)|0) . (35) 

These Green functions are deflned on a sphere of inflnite radius, at which the energy flow S{n) 



is measured |T^, |T3|, |n[. They depend only on the N unit vectors rTj, in addition to the two 
light-cone directions, p+ and p_, of the quarks, represented by the Wilson lines of Eq. (|2^) . In 
contrast to fri^), the Green functions Q are truely universal, depending only on the underlying 
dynamics rather than the choice of event shape. 

It is interesting to note that a stereographic projection maps each of the unit vectors n 
into a point on the complex plane p = {z,z), so that Q can be considered as a function on 
the 2-dimensional p— plane. In addition, Q is singular only when two or more nj's coincide, or 
equivalently the points pi approach each other on the 2-dimensional plane. These properties 
suggest that it possible to relate the Green functions ( ^5|) to (A^ + 2)— point correlation functions 
in some effective local two-dimensional fleld theory. 

Let us explore further the relation of the functions Q to the thrust, by considering the lowest 
moments of the thrust distribution Eq. ( [T^ ) 

it) = ^{e) + {tU , {t') - {tf = -^ Ub') - {Bf] + {t\, - {t)L . (36) 



g^ , . ^,PT, ^ , ^, Q 



^One can show that for the weight function w-{n) ~ 1 — cos6 the shape function /-(e) coincides with the 
inclusive (Hght-conc) structure function F{x) ||2^ in the hmit 1 — x = e/Q -^ 0. 
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Here, the moments of e are defined through (0) and (•••)px corresponds to the perturbative 
contribution. We observe that perturbative and nonperturbative effects contribute additively to 
the moments of the thrust distribution and that the 1/Q correction cancels in the variance of t. 
The leading nonperturbative corrections enter (^) through the average values of the powers of 
e in Eq. 



{e) = Ai = / dniwrini) G{ni) 
{e^) - (e)^ = -A2 = / dnidn2 wt(^i)w^t(^2) [ ^(^1,^2) - ^(^1)^(^2) ] • (37) 



We see that the mean value of e for the nonperturbative function /t(^), and hence the 1/Q power 
correction to (t) in Eq. (|36|) , is determined by the average value of Q{n), weighted by wxin). The 
variance of e, however, is related to the "connected" part of the Q{ni, 7^2) • 

Similar considerations apply to the heavy mass distribution. Using (^Tf) we calculate the 
lowest moments and obtain 



(P) = 7T7t(^) + (p)pt) 



1 
2Q 

1 



{p')-w = {p\.-{p)U + 



PT /m2 



4g2 



{e')-{ef+{{eR-eLf) {l + A dp'p 



daj 



PT\ 2 



dp' 



{31 



where the nonperturbative scales (e") = {{sr + El)"') {n = 1,2) are the same as for the thrust 
moments and a new scale {{er — SlY) is defined as an average with respect to the shape func- 
tion f{eL,£R)- We observe that the 1/Q correction to the mean value (p) is half that for (t), 
independent of the form of the shape functions. However, because of their different dependence 
on El and Er, starting from the second moment nonperturbative corrections to the thrust and 
heavy jet mass are not simply related, unless one chooses a particular ansatz for the underlying 
shape function f{EL,ER). 



The functions G defined in Eq. (35) measure the energy fiow in the final state, and it is not 



surprising therefore that the energy-energy correlation (EEC), 



1 

Q' 



EEC(x) = T^E^*^^- '^(cosx - cos%) , (39) 



«j 



has a simple interpretation in terms of these correlation functions. The leading 1/Q nonpertur- 
bative contribution to the EEC is associated with soft gluons propagating in the direction defined 
by the angle x, and it may thus be written as the expectation value of a single ^(n). Away from 
the jet directions, the soft gluon radiation that contributes to the EEC factorizes from the jets, 
and the EEC may be written as the usual expectation value for soft radiation from a product of 
eikonal lines, 

(EEC(x)),p = ^{0\W^{0) S{n) W^(0) |0) = l^(x) , (40) 

where x is the angle between n and the two-jet axis. We emphasize that this result holds only for 
X fixed and away the directions x = 0,7r, where Qsinx is a new, small parameter. Comparing 
Eq. ( ^OD and Eqs. ( P7| ) and ( |5B| ) we observe a relation between the leading 1/Q corrections to 
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EEC(x) and mean value of the thrust. This relation depends, however, on the angular form of 
the function G{x)- 

We can now relate our formalism to previous recent work on power corrections. Consider first 
the simplest ansatz for the ^'s, in which all higher-order connected parts of the Green functions 
Q are neglected, 

N 

gini,...nN) = l[G{ni). (41) 

According to Eqs. ( P^D and (^Tj), this is equivalent to replacing the shape functions fri^) and 
fi^Ly^R) by delta functions. 

Me) = 5(6 - Ai) , fisL, sr) =6(6l- ^Xi) 5 (^er - ^Ai) . (42) 

In this case, the variances of the nonperturbative distributions in e vanish, and substitution of 
( ^2]) into (p!9D yields a simple shift in the perturbative thrust distribution, as in Eq. (|l]). The 
same happens for the heavy mass distribution, but with half the shift in p. 

These event shape power corrections can be further related to the energy-energy correlation, 
if a form is given for the function g{n) itself. If we assume, for example, that the angular 
dependence of Q is identical to that of single-gluon emission diagrams, Q{n) ~ 1/ sin^ x, Eq. (|37| ) 
implies that 

Six) = -^V , (43) 

Arc sm x 

which relates the 1/Q power correction of the thrust to the EEC. In the model of 0, the 
parameter Ai is related to a nonperturbative, "universal" coupling constant a^^, although our 
analysis does not require such an interpretation. Comparison of the form Eq. ( ^Bj ) to data will 
help gauge the size of corrections associated with the variance and higher moments of the Green 
functions Q. 

There is evidence that corrections beyond Eq. (|^) are necessary. For the thrust distribution, 
the simple shift of Eq. (|l|) improves agreement with the data |15|, relative to resummed pertur- 



bation theory, for t > Aqcd/Q, but for smaller values of t ~ Aqcd/Q higher-order connected 
correlations evidently become important. The variance, for example, includes information on 
branching of final state gluons in all possible relative configurations, including opposite hemi- 



spheres [18|. Its contribution to (H2f) smears the 5— function distribution and gives rise to a 



"nonfactorized" correction to the shape function, Af{eL,eR) in Eq.(|TB|). As a simple ansatz for 

lese effec 

f^LSRY 



the shape functions that models these effects, one may consider the following expressions |TC 

a-l 




Me) = [-) exp^-^J^. (44) 

with a and A free parameters. Their values can be fixed using (|37D as 



A^^^^^ = Ai , A'a = \l-X,. (45) 

T{a) 

A fit of this kind for fx, using data at a single value of Q - the Z mass - induces a good description 
of the data for the thrust distribution for energies 14 GeV < Q < 161 GeV, over the full range 
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< t < 1/3 as shown in Fig. 1 [0, ^. Further analysis of event shape energy dependence 
should make it possible to estimate the underlying energy flow functions Q{ni). 
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Figure 1: The comparison of the data with the QCD prediction for the thrust distribution at 
different energies (from bottom to top): Q/GeV = 14,22,35,44,55,91,133,161, based on the 



shape function. The detailed description of the plot can be found in [T^ . 



5 Summary 

In this paper we have studied the power corrections to the differential thrust, t = 1 — T, and 
heavy mass, p, distributions in e"'"e~ annihilation close to the two-jet limit. In addressing this 
problem, we did not aim to justify a particular QCD-inspired phenomenological model, but rather 
to formulate a framework with which to study the relationship between perturbative and nonper- 
turbative effects in high-energy final states. We have seen that, despite the fact that the thrust 
and heavy jet mass are not inclusive quantities, the leading nonperturbative corrections to their 
differential distributions can be factorized into the perturbative and nonperturbative functions, 
in much the same way as for inclusive cross sections. We identified nonperturbative infrared 
shape functions that organize all leading power corrections, l/(t(5)" and l/{pQ)"'. Although not 
universal themselves, these shape functions can be derived from universal matrix elements that 
describe energy fiow. We anticipate that it will be possible to extend these considerations to a 
wide class of infrared safe event shapes and hard-scattering processes. 
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